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Abstract 

A nonlinear sigma model is derived for the time development of a Bose-Einstein condensate 
composed of fermionic atoms. Spontaneous symmetry breaking of a Sp(2) symmetry in a co- 
herent state path integral with anticommuting fields yields Goldstone bosons in a Sp(2)\U(2) 
coset space. After a Hubbard-Stratonovich transformation from the anticommuting fields to 
a local self-energy matrix with anomalous terms, the assumed short-ranged attractive inter- 
action reduces this symmetry to a SO(4)\U(2) coset space with only one complex Goldstone 
field for the singlett pairs of fermions. This bosonic field for the anomalous term of fermions 
is separated in a gradient expansion from the density terms. The U (2) invariant density terms 
are considered as a background field or unchanged interacting Fermi sea in the spontaneous 
symmetry breaking of the SO (4) invariant action and appear as coefficients of correlation 
functions in the nonlinear sigma model for the Goldstone boson. The time development of the 
condensate composed of fermionic atoms results in a modified Sine-Gordon equation. 
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1 Introduction 

Experiments of Bose-Einstein condensation (BEC) with bosonic constituents have been realized 
under various conditions. In many cases the Gross-Pitaevskii (GP) equation with a bosonic field 
tpg(t) as the order parameter and wavefunction can be applied 

lh ^t^L = _|_v2^(i) + u (x,t) &r(t)+2y)|^(t)| 2 Vj#_^ ips{t), (1) 
Ot 2m 

x' 

where u(x, t) refers to a time dependent external potential, including the trap potential, and Vj^_gi 
is a short-ranged interaction ^1 [5] . Transfcring the GP equation to the case with fermionic 
atoms, a coherent field equation with anticommuting numbers can be introduced where the classical 
field ips{t) is replaced by a Grassmann-valued field xs.s(t) with spin s =|, J, [H] 

lh ?XjM = _|lv 2 X*,.(*) +«(*,«) **,.(*)+ (2) 
Ot 2m 

+ 2 X*S',a'(t) %-2] XS',s>(t) XsAt)- 
x' ,s ; 

As its bosonic counterpart, the Grassmann-valued equation is integrable for a contact interaction 
and possesses a set of infinite independent integrals of motion jlj . This has been demonstrated by 
the method of Lax-pairs and r-matrix methods. 
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In this paper an effective nonlinear sigma model of bosonic fields for a condensate composed 
of fermionic atoms is derived for the following quantum Hamiltonian of Fermi operators ips,s 
corresponding to the classical equation J2J with anticommuting fields Xx,s(t) 2 

H^+^t) = V£ jB / Vfc,. (3) 

x,s;x f ,s' 

+ E ^ta V \S'-S\ ^x,s i>3>.s> 

x,s;x' ,s' 

X 

Hs',s'-s,s{t) = S s / jS Sg/ > g(^—+v(x,t)\ (4) 

v(x,t) — u{x, t) — fiQ . (5) 

The external potential u(x,t) is shifted by the transformation exp{— i/h ■ fj,o t} ip£^ s with the 
chemical potential /iq. This shift is performed because the time derivative is considered as a 
perturbation in a gradient expansion so that rapid oscillations of the fields do not appear. Since 
we also include a time dependence in the external potential u(x,t), the chemical potential /j,q 
is in general not an equlibrium value and can be extended with an appropriate adiabatic time 
dependence fi(t). However, we assume that the interacting Fermi sea is not strongly perturbed 
by the time dependence of u(x,t). The interaction V\g>-g\ < is attractive and has to be short- 
ranged in order to obtain local sigma matrices for the self-energies. Furthermore, the densities of 
the interacting Fermi sea are regarded as given background fields which are considered as given 
coefficients for the gradients of the bosonic nonlinear sigma model. The spatial gradient expansion 
for the nonlinear sigma model is combined with a kind of Thomas-Fermi approximation 2. where 
the derivative terms of the kinetic energy are taken into account as a perturbation [S]-j7]. In the 
remainder we investigate and assume a BCS like condensation phenomenon of the atoms, derived 
from spontaneous symmetry breaking with the source field jg(t), and exclude the formation of single 
bosons from bound pairs of atoms because the attractive potential is taken sufficiently short-ranged. 
In the case of a box potential with depth Vq < and range r , this means that m \Vq\ r^jh 2 < tt 2 /8 
has to be sufficiently small 8 . This case with short-ranged attractive interaction is different from 
the formation of excitons and biexcitons of the long-ranged Coulomb potential in semiconductors 
where the Hamiltonian for the strongly bound electrons in semiconductors can be reduced to 
purely bosonic operators jOj. A review of the nonlinear sigma model in superconducting systems 
with delta-function correlated disordered potentials and the replica-trick can be found in Refs. 

maim. 

A coherent state path integral with Grassmann fields Xs,s{t P ) = Xy(tp) (H = s }) 3 is 

used on a nonequilibrium time contour J c dt p . . . = dt + . . . + J^™ di_ ... to express the time 
development of the system with Hamiltonian © 



Z[J] = J d\xy(t P )} (6) 
6XP { ~ ^ Jc dtp 5 X * W ^ Sff '^( ~ ~ * £p ) + Hff '^ tp) Xff ^} 



!/■■!/ 



Cxp { ~ k j dt P^2 x *y'^ x *v^ V \2'-A Xa(tp) Xy'(t P ) 
c v;v' 



2 The spatial sum y]_ ... is dimensionless and is scaled with the system volume so that ... is equivalent to 
j) , A d x/L d .... 

3 In the following {/-vectors refer to the combined spatial vector x and spin variable s =f , J, as abbreviation. 
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x expj - ~ J^dt p ^2jt(t p ) (xs,r(t P ) Xx,i(t P ) - Xxditp) X?,f(*p)) +h.c.| 

Xa(ip) = X?, s (ip) (7) 
Wfe) = =( 2'!(*p) ) (8) 

A source field jg(t p ) is applied to generate bosonic pairs J6?,t(*p) Xs,J.(*p) ou t of the vacuum 
state for spontaneous symmetry breaking of the U(l) invariant one particle part 7i and the U(l) 
invariant interaction [201 E]- In order to generate observables, the source term J$>-$(tp} has 
to be introduced where the fields xs,i(t P ) and Xx,i(tp) have to be combined to the four component 
vector rjg_ s {t p ) JHJ so that pair condensate terms X3',s'if,p\ ) X£,s{tpi) can be obtained by simple 
differentiation with respect to the matrix 3y'-$(t v 2 }\ tp 1 ^) |22| . A nonhermitian infinitesimal part 
— i e p — —i (±e), (p = ±) on the time contour has to be included for the analytic and convergence 
properties of Green functions, derived from the coherent state path integral Z[J] (JHJ. 

Apart from the commutation relations the field equation, Hamiltonian and coherent state path 
integral are formally similar in terms of Grassmann numbers to the bosonic case |16| . Therefore, 
variations of the action in JJjJ and other approximations can be performed, however, compared to 
the condensation of single bosonic constituents, only a small fraction of fermions can condense near 
the Fermi energy pQ . This means in terms of a density matrix formulation that expectation values of 
densities (Xx',s'{tp) X*s a (*p)) are considerably larger than the pair condensate (Xx'.s'{tp) Xx,s(t p )}. 
It is the aim of this paper to extract the various densities and pair condensate functions from the 
coherent state path integral and to derive effective equations for the pair condensate composed of 
fermionic atoms, in analogy to the GP-equation for bosons or fermions The nonlocal order 
parameter ^y'-y(t p ) has a matrix form 

%-';y(t P ) = ( g'^j ) ® ( Xi,M XxAh)) (9) 
f (Xx',s'(tp) X*S. s {t p )) (Xx',s'(tp) XxAtp)) \ 

V b&,A t p)**M) -(xzAtp) x* s ,,Atp)) /, 

where a doubling of the spin space has to be considered because of the source field j~(t p ) which 
causes the spontaneous symmetry breaking. The order parameter is invariant under U(2) 
transformations in spin space which does not alter the block structure into densities and pair 
condesates. The form of the order parameter also allows a global hyperbolic symmetry which 
combines densities and pair condensates. A complete symmetry group of the path integral is 
spontaneously broken by the subgroup U(2) for the invariance of the densities and the source 
term. This symmetry breaking leads to a nonlinear sigma model after a gradient expansion for the 
anomalous terms. The various steps for obtaining the nonlinear sigma model are briefly listed : 

• coherent state path integral 

• transformation of the quartic interaction to densities, anomalous terms and the order param- 
eter 

• Hubbard- Stratonovich transformation from the fields to the self-energy and integration over 
the remaining bilinear anticommuting fields |2.'il 

• the short-ranged attractive interaction reduces the Sp(2) symmetry of the path integral with 
a spatially nonlocal self-energy to a spatially local self-energy matrix with 50(4) symmetry 



x exp 



i 

2h 
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• separation of the self-energy into densities and non-diagonal terms on a coset space according 
to spontaneous breaking of the orthogonal symmetry SO(4)\U(2) and determination of the 
measure 

• separation of the coherent state path integral into block diagonal U(2) invariant density 
matrices and anomalous terms including a gradient expansion 

2 Hubbard- Stratonovich transformation and self-energy 

The quartic interaction with the short-ranged two body potential has to be transformed to a relation 
with an order parameter similar to &y'-,y(t p ) @. The antisymmetric fields Xx,s{tp), s (t p ) and 
Xx',s'(tp), X%i s '(tp) can De combined in the following even matrices rg'.g(t p ) = rg/ s >.g tS (t p ) and 
Py';y(t p ) — ps> .s' :£,s(t p ) where r and p are hermitian and antisymmetric, respectively 

r y';v(tp) = Xy'itp) X*y(tp) (10) 
r*tf;y{tp) = Xy(tp)x*y'(t P ) = r ff , ff/ (t p )^r + (t p )^r{t p ) (11) 

Py';y(tp) = Xy> ih) Xyitp) = ~Xy{tp) Xy> ih) = ~ Py:y> (*p) ~> p(tp) = ~ p T '(*p) (12) 
Pk'Sv) = X*y( t p)X*y'(tp)=P^(t p )=pi ff ,(t P ) ■ (13) 

In terms of the even matrices r and p, the quartic interaction can be written in the form of an 
order parameter R~ b s .gr s i(t p ) as in with a doubling of the dimension of spin space where the 
superscripts a, 6=1,2 refer to the doubling and s, s' =|, | to the spins 

^ X*S',s'{ i 'p)X*x,s{ t p)Xx,s{ t p)XS',a'{tp)V\g>-g\= (14) 



x,s\x' ,S 



4 (Xff'itp) Xy'itp) ~ Xy'itp) Xff'itp)) (x*yitp) Xyitp) ~ Xyitp) X|"(*p)) V \S'~S\ 
y;y' 

~ 4 Tr Rg b S ;x' ,s r itp) ( 2 J ^iP,s';x,s(tp) ( 2 _y 2 



Rfs-^Atp) - ( g;^) ) ® ( xxv(^)) as) 

Tx,s;x' ,s' itp) Px,s;x' ,s' \tp) 
Px,s:x' .s'itp) ~ r x,s;x',s'itp) 

Obviously, the quartic interaction can be expressed with the nonlocal order parameter Rg b s .gi s iit p ) 
so that a global hyperbolic symmetry results with the diagonal matrix k = diag(l, 1, —1, —1) and 
the matrix T between densities and pair condensates (T + kT = n), apart from a [/(2) invariance 
in spin space 

R -> TRT+ (16) 



T = ( Vi r" . /T~7~77T ) * = * T (17) 



Vi + t+t t+ 
t Vi + tt H 

T+ ( 12 - la ) T =( 12 -ij ■ < 18 > 



The 2x2 matrix t ss > in spin space of the 4x4 matrix T s a s h , has to be complex symmetric and 
therefore contains 6 real parameters (see appendix 0) . The hyperbolic symmetry with matrix 
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T and the U(2) invariance in spin space is equivalent to a symplectic symmetry group Sp(N/2) 
and not a unitary group as U(N/2, N/2) because of the number of independent parameters which 
equals ten in the considered case. The symplectic symmetry becomes obvious after an exchange 
of the first, second with the third, fourth row of the diagonal matrix k and reordering the matrix 
T + to its transpose T T in relation 1)18(1 (see appendix 0) . Since the matrix T has dimensions 
4x4 and the number of independent parameters for Sp(N/2) is ^N(N + 1), a Sp(2) invariance is 
obtained for relation 1(14(1 . Using the identity for the Hubbard-Stratonovich transformation with 
the self-energy matrix SS s ,(t p ) consisting of commuting elements only 



^x,s:x',s' ftp) 



Sx,s;x' ,s' ftp) ®x,s;x',s' ftp) 
<J x,s;x',s' ftp) ~ S x,s;x',s' ftp) 



(19) 



the quartic interaction term with V[j/_j| can be expressed as a quadratic term of the self-energy 
and a bilinear product of anticommuting fields rjs, s ftp) = (xs,sft P ), Xg s ftp)) T 



exp j - | J dt p J2x*y<(t P ) Xffftp) V\x'~x\ Xyftp) Xy'ftp)\ = J d W d[a] 
c y;y' 



" X|, { ^I' lll ^.V7-J-. 

my' 

k = diag(l, 1,-1,-1) 



Xy'{tp) 



(20) 



Substitution of the interaction with the Hubbard-Stratonovich transformation yields the following 
coherent state path integral where a doubling of the one particle terms has to be performed because 
of the source terms so that a bilinear product of anticommuting fields is obtained 



Z[J}= I d[s}d[a] cxp\-^- [ dt P Y, ' 



Tr 

v\t' — t\ s.s'a.b 



S K £ K 



(21) 



d[x] exp 



2h 



dt p r] + 



-J + -H T 



J + S 



3 88 f l^p) 



_0 j(x,t p ) 
-j(x,t p ) 
d 



P 



H ss >(t p ) = S ss > ( - l ^-Q t is p + — + v(x,t) 



(22) 
(23) 



After integration over the anticommuting variables, the path integral only contains the self-energy 
£ with the hermitian matrix s for the density terms and the antisymmetric matrix a for the 
anomalous terms 



Z[J] 



d[s] d[cr] exp I [ dt P Y^ 



. Tr 

- -., *\3'— S\ s,s' ,a,b 



(24) 



\ 



det 



nt P ) -j{tp) 

-j + (t P ) -n T (t p ) 



J 



s{t p ) -cr{t p ) 
-a+(t p ) -s T (t p ) 
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The generating function Z\J\ l)24|) is invariant under a symplcctic symmetry group Sp(2) which 
has its cause in the anticommuting properties of the fields xs,s{t P ) and the doubling of spin space 
(see appendix IXjl . This symplectic invariance is spontaneously broken by the U(2) invariance of the 
matrix s and — s T and the source term j(t p ) l122[| so that three complex or six real Goldstone fields 
result on the coset space Sp(2)\U(2) because the corresponding dimensions of the Lie algebras for 
Sp(2) and U(2) are ten and four, respectively. The densities s and — s T represent the self-energy 
of the interacting Fermi sea and can be regarded as a vacuum state or background field on which 
the subgroup U (2) invariantly acts so that the symmetry Sp{2) of the complete Lagrangian is 
spontaneously broken to three complex Goldstone fields. 

However, if the trap potential in v(x,t) can be restricted to a typical distance clq and if this 
distance ao is considerably larger than the range ro ~ \x' — x\ of the two body potential V\x'~x\i 
there are strong oscillations in the quadratic term with the nonlocal self-energy 'Efs-x' s , {t p ) of (|24() 
because f/V\g'-2\ tends to infinity as the short ranged potential V|^'_j| approaches zero. Therefore, 
the spatially local parts E^, (x, t p ) of the self-energy are only retained in the path integral (|2"4f 

— — — — * oo for \a? — x \ > r (25) 

*\x' —x\ 

^Z b s;x' ,s' (tp) ~ * ^ss'(^!*p) &S',$ ■ (26) 

This can be accomplished by integration over the matrix elements E^.-, ,(t p ) with \x — 5?\ > ro 
in l|24l) and eliminates the nonlocal parts in the self-energy which cause the strong oscillation 
on the nonequilibrium time contour in (12411 . The diagonal parts crg-i-sP t(t p ), erg i;x',l(tp) of spin 
space in the matrix o~(t p ) tend to zero because of the antisymmetry of the pair condensate in the 
spatial part. The vanishing of the diagonal elements of spin space in a(t p ), due to the assumed 
short-ranged and spin independent interaction V\g^gi \ , corresponds to the observation that there 
is usually no triplett pairing of fermions in the condensate. Consequently, only one complex 
Goldstone field for the singlett mode remains whereas the other two complex fields, which result 
from spontaneous symmetry breaking Sp(2)\U (2) in the path integral with nonlocal self-energy, are 
suppressed because of the short-ranged interaction. 4 A Sp(2)\U(2) coset space for the Goldstone 
bosons would result if the interaction V\s'-s\ was constant for any pair of spatial vectors x' , x so 
that every atom would interact with equal weight with all other atoms independent of distance. 

After a shift of the self-energy matrix by the spontaneous symmetry breaking term, the coherent 
state path integral Z[J] (|24ll is transformed to a local self-energy T,^,(x,t p ) consisting of the 
hermitian local matrix s and the antisymmetric local matrix a in spin space and the two body 
potential restricted to a finite typical zero distance value Vq < 5 

Z[J] = J d[s] d[<r] (27) 



exp 



exp \ / dt p y~* -jr Tr 

4ft J C -V s,s',a,b 




The 4x4 local self-energy matrix in (|27|l consists of the hermitian U(2) invariant density term 
s(t p ) with four real parameters and the antisymmetric 2x2 anomalous term with one complex 



field so that the matrix (x, t p ) contains six real fields. The number of independent parameters 



4 In the case of spin dependent forces or other interactions Vjg_j/| which have their maximum for \x — x'\ ^ 0, but 
vanishing zero distance interaction Vq, other complex Goldstone fields have to be chosen. These cases are excluded 
in the present paper. 

5 The variable M in 1271 is a normalization factor J\f = (L / Ax) d ■ (1 / At) because a determinant without integration 
measure is considered in the Tr In term. 



(> 



and the dimension of Jl^,(x, t p ) indicate a S'0(4) symmetry where the antisymmetry of S'0(4) 
generators becomes obvious after exchange of the first, second with the third, fourth row of the 
4x4 matrices in (|27H . The following local parametrization (|31I33|) of the self-energy with only one 
complex field <p(x,t p ) as Goldstone boson can be chosen in the coset space SO(A)\U(2) in order 
to separate the anomalous term from the unchanged interacting Fermi sea with density matrix s, 



in the spontaneous symmetry breaking 



^ss' i X i tp) ~ 

S*/ S (x, t p ) = 
<7 TT (£, t p ) = 



(x, tp) 

0~s's 





r + 

ss 

Sss' (*£ ; tp) 

° 'u(^) f p) 



Q ss' \<E<, tp) 
i (x, tp) 



(x, tp) 



i (x, tp 



^ss' ( x i tp) 



T[JU ' ,] ( - SD( o' ~sl(x,t p ) ) 



T 
t 



t+ 



VT+t+t 

t y/l + tt+ / 

X, t p ) I2 W = W r + I W- 





' u 


w \ 


) •■>-{ 




V j 









u, v, w r , Wi € R 



(28) 

(29) 
(30) 

(31) 

(32) 
(33) 



The coherent state path integral can be transformed with the chosen parametrization (|31I33|) of £ 
and the change of integration measure w?(x,t p )/4 to the form 



Z[J] = j d[u] d[v] d[w r ] d[ Wi \ d[</>] ( n 

■ {S,t p } 



wf(x,t p ) 



(34) 



CXP {-4k/ c ^5^( 2 ^^- 4 ^- 16 
exp \ ~ j dt p y\M Tr In 

* JC - s,s',a,b 



■ j) VT+W- 



J + T 



sd 



D 



T+ 



Since only a small fraction of the Fermi sea condenses, classical equations for the 2x2 density 
matrix Sd(x, t p ) can be obtained by variation of the action in l)34|l with respect to u, v, w r and 
Wi (1321331) where the matrix T is set to the unit matrix and the integration measure wf{x,t p ) is 
included in the variation 6 



variation Su(x, t p ) : 

-J-u + (H + s° D )- 1 (x,tp)=0 

Vq 

variation 8v(x, t p ) : 

v° + (H + s° D )- 1 (x,tp) = 



V 

variation 5w r (x,t p ) : 

~ w° r + (H + s° D )^(x, t p ) + (H + s° D )j{(x, t p ) = 
variation Swi(x,t p ) : 

11 ' kn + s%)rf@, tp) -(H + s° D )ll(S, tp) 



(35) 
(36) 
(37) 
(38) 



V 



w° + 1 



2 

to? 



To =° 



The classical equations (I35I38|I for the resulting matrix s%(x, t p ) can be simplified with the Thomas- 
Fermi approximation where the kinetic energy can be neglected because of the large atom masses 

6 The fields u(x,t p ), v(x,t p ), w r (x,t p ), Wi(x,t p ) in S£>(x,t p ) have to be scaled by the factor (At/h) ■ (Ax/L) d 
to dimensionless quantities for the variation. 
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and an assumed homogeneous self-energy of the bulk Fermi sea. This gives algebraic equations 
for Sjj(x, t p ) which can be applied in the correlation functions of the nonlinear sigma model with 
matrix T (|32[1 following in section [3] Using the parametrization into block diagonal densities 
for the Fermi sea and anomalous terms, the determinant in (|34fl has to be expanded with respect 
to the gradients contained in H [S]-[7|. 



3 Separation into densities and anomalous terms with a gra- 
dient expansion 



Applying the chosen parametrization 1)31133(1 for S, we can insert the term k T k T + k into the Tr In 
term of l|34l) without a modification of the coherent state path integral because the determinant of 
k in combined spin and hyperbolic space equals unity 



Ox = [ dt p VA/" Tr In 

JC - s,s',a,b 

x 

= [ dtpS^JV Tr In 

J C - s,s',a,b 



n 



-n 1 




(39) 



H + sn 



n 1 + s 



T k J T~ 



T 



H 



n 1 



T -i_ 



n 



n 1 



SH a s b s , 

The gradient expansion of STif^i gives the following operator S7i" b , — S ss i 8H ab 



T 



H 



n 1 



T- 



H 



TL 



(40) 



(T k T- 1 - k) {-Ep) + T k (-EpT- 1 ) 



2m 



E,„ = ih 



—dp (dpT) T- 1 
2m M y ^ ' 



0_ 

dt„ 



(41) 



where one has to distinguish between the operators dp, E p and the derivatives (dpT),(E p T) = 
ih{dT/dt p ) of the matrices T, T~\ e.g. T dp T" 1 = T UdpT^ 1 ) + T^ 1 dp . 

Expanding the Trln term 0\ up to second order in SH^ b ,, we obtain the expression (|42[1 with 
the Green function G a ss , of the block diagonal density matrix sd in spin space in symbolic form 
(spatial and time coordinates are omitted for brevity). In the following the field Wi in the matrix 
sd Q32I33JI for the Green function G® s , has to be separated from the expansion because it couples 
to the orignal U(l) symmetry violating source term j I(34I6|I 



(42) 



d = 2 / dipVA/'tr ln(n + s D 
Jc 3 v 

dtp^M Tr \G a ss , (T kJT- 1 + SH) 



s,s ,a,b 



- / dt^Af Tr 



s.s ,a,b 



Gt, SH ab G b s , s SH ba 
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£ J n _ s.s ,a,b L i 122 

x 

dt p ^TAT Tr \G ai , (T n J T~ 1 ) a } a ? G a ? 5H a2t 

X 

The hermitian self-energy SD-ss'(x,t p ) without the field Wi is only contained in the Green function 
G* a , (x, t p ) x', tp) on which the derivative operators in 8H ab act. A spatial and time diagonal Green 
function g ss i (x, t p ) for G" s , can be considered in a Thomas- Fermi approximation for large atom 
masses and a nearly homogeneous system where the kinetic energy and time derivative can be 
added as a perturbation 

G a s { s 7 1/2) (2, t p ; Sf, t'p) = (x, s; t p \(H + Sd )( t >- V , (43) 
= I — ie p + v(x, tp) + sd{x, tp) ) Sg^i S(t p — t') + derivative terms . 

V / ss' 

v ' 

The first order term of STC^, in 0\ (second term in I42JI vanishes in the Thomas- Fermi approxima- 
tion with G a ss , replaced by the spatial and time diagonal function g ss i(x,t p ) l|43[l . The anomalous 
terms with the matrix T can be reduced to a 2 x 2 matrix because the spin space is restricted to 
the block diagonal densities. 

Introducing the following averages of the block diagonal densities with matrix SD-ss'{x,t p ) 



d[s D ] n det(n + SD ) 



{«,*„} 



(44) 



Wo Cu(x,t p ) 



c t (x, tp) 



i -* \ 

W c ^ x ^p> 



hv 



c(x, tp) 



— — C fi (X, tp) 

nvo 



hv 



Cf_ti> (% 'j tp) 



(E p g as >(x,t p )) (E p g s , s (x',tp)) 
(-E p g S s'(x, t p )) g s ' s (x,t p ) 
(d M g as '(S,t p )) (-E p g s > s (x,t p )) 

Qss' tp) Qs's (^5 tp 

g ss '{x,t p ) (dpg S ' S (x,t p )) 
(d,tg ss >{x,tp)) {d v g a >s(x,t v )) 



(45) 
(46) 
(47) 
(48) 
(49) 
(50) 



the coherent state path integral can be separated into an action 5o[s_d] following from (|44J) and an 
action S[T, T -1 ; {cij}] for the anomalous terms 



Z[J] 



d[ SD ] n 



{x,t p } 



cxp I - i 5 [s D ]| 



(51) 



J d[0] exp { - i S[T, T~ l ) {Cij}] } exp { - i Sjj [T, s D ;J, j] } 



The action S[T,T ; {cy }] is given by the following relation up to second order in the gradients 
(9 M , E p where partial spatial integrations have been performed. The parameter functions Cij(x,t p ) 
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(|45I50J) contain the properties of the densities with the matrix as a background field 

T" 1 = kTk k = diag(l,-l) 

{VT- 1 ) = -(EpT- 1 ) = ih T' 1 (8 tp T) T- 1 



S[T,T- 1 ;{c ij }] = ±-J c dt p J2 



(52) 
(53) 

(54) 



c(x,t p ) tr {Tk{T>T )) - CliV (x,t p ) — tr (d^T) (d v T ) 



,2x2 



r i 2h r 

+ c tt (x,t p ) tr (TkT^-k) 2 + Clx {x,t p ) tr (d^T)K(VT^ 1 ) 

a.b L -I Tfl a.b L 

+ (2 c, lt (x, t p ) - OpCtix, t p j) — tr \[T, k] (dpT) 

V / Tfl a,b L 

+ 2c t (x,t p ) tr [(T _1 — T) (VT) 



It can be verified with an expansion of the Goldstone field <fi(x, t p ) of the matrix T in sinh-amplitude 
and phase term <fi(x,t p ) — smh(ip(x,t p )) exp{i a(x,t p )} (ip(x,t p ) > 0) that the derived action 
(|54|) is composed of a massless Goldstone field a(x, t p ) and a real massive field smh((p(x,t p )). 
The parameter functions Cij(x,t p ) can be regarded as generalized mass and kinetic terms of a 
spontaneously broken 4 field theory [21]. They can be calculated from the classical equations 
()rt5lrt8JI for s° D (x,t p ) where the Thomas- Fermi approximation of large atom mass can be used for 
simplicity. In terms of the sinh-amplitude and phase term, the action S[T, T ; {c^ }] l|54|l is similar 
to that of generalized Sine-Gordon equations 



% ' a;{Cy]1 = i/ c d '^ 



(55) 



2h 2 sinh 2 (<^) 



c(x,t p ) cosh(2<p) (d tp a) + c tiu (x,t p )[ — ) (d^a) (d u a) 



2fi z 



2h 



c(x, t p ) (d tp tp) 2 + Cp V {x, tp) 



(d^ip) (d v tp) 



- 8 c tt (x,tp) sinh 2 (</?) 



+ c M (f,t p ) sinh(2^) (d tp a) (dptp) - (d^a) (d tp <p) 



Aih 1 



{2 c^ t {x, t p ) - d^c t (x, t p )^j sinh 2 (tp) (d^a) + A1T1 c t (x, t p ) swh(2tp) (d tp tp)j 

A first order variation of the action i|55|l with respect to the fields on the time contour gives a classi- 
cal equation for the time development of the Goldstone field 4>(x, t) = sinh(^(.x, t))x exp{i a(x, t)} 
with the densities SD;ss'(x,t p ) as background fields in the parameter functions Cij(x,t) 



tp(x,t) = ~(tp(x,t+) + <p(x,t^ 
a(x,t) = — (a(x, t + ) + a(x, t_)^ 
Cij(x,t) = -(cij (x,t+) + Cij (x,t-jj 

d t (c sinh 2 (^) cosh(2^) (d t afj + ( — J d^{c pv sinh 2 (^) (d„a) 



(56) 
(57) 
(58) 

(59) 
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+ ^ smh(2(p) (^(dtCft) (d p ip) ~ (<^ c m) (9^)) - —9^(2^ - (c^c*)) smh 2 ^)^ = 

9 t (c (d&)) + (£\ {djip)) = (60) 

sinh(2<^) x |c(cosh(2^) - \)(d t af + ^L(d^a)(d v a) - ^c tt - -(efe) + 

The equations 1)59160(1 determine the time development of the Goldstone field for the condensate 
composed of fermionic atoms and are analogous to the Gross-Pitaevskii equation for bosonic con- 
stituents. 



4 Summary and discussion 

Since the parameters Cij(x,t p ), following from the densities in the background, change slowly in 
time and spatial coordinates, the pair condensate is determined by the lowest order terms of the 
nonlinear sigma model with the matrix T 



T= (V^+W J_) (61) 

This separation into background properties with the parameters cy and anomalous terms of the 
nonlinear sigma model corresponds to the observation that only a small fraction of the atoms 
condense in the Fermi sea. In the case of computations one therefore can introduce correlation 
functions of the densities, as e.g. the current-current correlation function c pv (x,t p ), etc. which can 
be calculated with the self-energy matrix s° D (x,t p ) 1(35138(1 of the first order variation of the action 
in (|34l) . Taking derivatives of the action Sjj[T, sd', J, j] with respect to J, the appropriate ob- 
servables can be obtained, as e.g. the anomalous term -jf(Xx,i(t p ) Xx',i(t p )) of two anticommuting 
fields is represented by a relation with the bosonic matrix T and the Green function G^ s , of the 
densities 

-^jf(XsAtp)Xs>M))= ( 62 ) 

2 

= ^^(T 1 ) 2a (x',tp) G a ^{x\ t' p ; x, t p ) — G"^ (x, t p ; x', t' p ) T al (x 7 t p ) . 

a=l 

Statements about binding energies of the pairs (Xx,i(t p ) Xx,i(t p )) can be calculated by various 
saddle point considerations from the coherent state path integral ((24(1 25, 26 . However, the exact 
expression of the saddle point is not needed for the derivation of the nonlinear sigma model with 
actions ((54(1 and l(55|) for the time development of a condensate composed of fermionic atoms. The 
equations 1(591) and 160() simplify considerably for a translation invariant system with a spatially 
constant and time independent external potential u(x, t) JHJ) so that only terms with the constant 
coefficients c ((48(1 and ((50(1 remain. In this case a simpler form of the Sine-Gordon equation 
results for the condensate of the Goldstone field 4> in the matrix T 1(61(1 . 

In experiments fermionic 6 Li was cooled to degeneracy by 23 Na or by mixing two different 
states in an optical trap (271 128j and fermionic 40 K was cooled to T/Tp = 0.3, (Ep — ksTp) by 
sympathetic cooling with 87 Rb Apart from degeneracy a BCS transition, for which the time 
development of the condensate is described by Eqs. ((44(1 to 1(60(1 . is suggested in Ref. [23j. For 
this a degenerate Fermi gas with attractive interaction between the fermions must be prepared. A 
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possible realization can be obtained in fermionic 6 Li 30 . In experiments atoms can be prepared 
in two different hyperfine states and a Feshbach resonance is used to tune the interaction due to 
elastic collisions to the desired value [3T]. In Refs. IHS-EHi BCS transition temperatures Tbcs 
are predicted in the range from Tbcs/Tf = 0.025 to Tbcs/Tf = 0.4, but it is an open question if 
these conditions can be reached experimentally. 

Close to a Feshbach resonance, the additional phenomenon of a BCS to BEC crossover is 
predicted for atomic Fermi gases |35) . A strong attractive interaction can arise between fermions 
mediated by bosonic quasi-molecules associated with a Feshbach resonance. In this case it is 
stated that one has to be careful in applying pure BCS theory to a Fermi gas when the pairing 
interaction is very strong |35|. due to fluctuations in the two-particle Cooper channel. One can 
also try to extend the experiments in optical lattices with bosonic constituents to fermionic atoms 
|36| . However, optical lattices involve the additional length scale of the periodic potential so that 
one has to examine whether a gradient expansion up to second order is sufficient for describing the 
atoms interacting on one lattice site. In this case one has to take into account the band structure 
of the optical lattice in the expansion of the Tr In- term of relation l|39|) . 

In the present paper we have performed symmetry considerations on which the derivation of 
the nonlinear sigma model with the Goldstone field <f> in the coset matrix T Ij61|l is based. Since the 
Sine-Gordon equation allows nontrivial solitonic solutions in 1+1 or two spatial dimensions, one can 
also expect solitons in condensates composed of fermions and study their temporal evolution. Apart 
from numerical computations, the derived effective Gross-Pitaevskii equation for fermions gives rise 
to investigations for Backlund-transformations in 1+1 or two spatial dimensions so that nontrivial 
solutions can be obtained as for the nonlinear Schrodinger equation of bosonic condensates. 
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A Symplectic symmetry of the generating function 



In the following we abbreviate spatial and time coordinates x, t p of the fields xs,s{tp), X*s s (^p) with 
the indices i, j in the global transformations acting in the doubled spin space. The symmetries 
result from the spin independent one particle Hamiltonian Ti and the spin independent interac- 
tion. Considering the Hamiltonian 7ijj as a spin independent matrix, we can rewrite the bilinear 
Hamiltonian form with the antisymmetry of the fields Xi.s, X* s as the relation 



Xi, a Xj,s = \ (xl,s H-H Xj, s ~ Xj,s 7-lij Xts) (63) 



i ( xts \ T ( % 




2 V Xi,« J \ 

After a transformation with a matrix M , one has to obtain again a separation into fields x'i s an d 
their complex conjugates x*' s m the first, second and third, fourth row. Therefore, the global 
transformations with the 4x4 matrix M consist of two 2x2 block matrices A, B in spin space 
and their complex conjugates 



Xi, s > \ ( A s , s B s , s \ / Xi,a 
X£ s , J \ B*, s A* s , s ) \ xls 



(64) 



M ty 



Because of relation l|63|) the matrix M has to be invariant under the following unitary hyperbolic 
transformation 

M+ kM = k k = ^ 12 _ 1 ^ (65) 

or the equivalent symplectic form which results from interchanging in k the first and second row 
with the third and fourth row, respectively, and from reordering the matrix M + to its transpose 
M T 

fT „ n/ r _ „ „ _ ( -12 



M 1 gM = g g = ^ ^ j (66) 

The bilinear form x* tXi,T + Xi \Xi,l i s invariant under these transformations with the matrix M, 
both for the unitary hyperbolic and symplectic kind of transformation 

2 {Xi,tXi,t + X*,iXa) = vt K Vi = Vl 9 Vi (67) 

/+ / IT 1 

= 2(x;: T x- T +xuxu) • 

The complex 2x2 matrices A, B consist of sixteen real parameters which are restricted by the 
equations (|65|l or 166(1 . yielding six real conditional relations. Therefore, ten independent param- 
eters remain in the matrices A, B. This indicates a symplectic symmetry group Sp(N/2) with 
7}N(N + 1) parameters for N — 4 of the 4x4 global transformation matrix M. The matrix M 
can be expressed with the global U(2) subgroup in spin space and the coset space Sp(2)\U(2) 
which is composed of the complex symmetric matrices t ss i or t' ss , with six and ten real parameters, 
respectively 



VI + 1'+ t> t'+ 



t' Vi + 1< t'+ 

U T tU t = t T t' = t' T . (69) 
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The matrices A, B and their complex conjugates are therefore given by the following equations 
which contain ten real independent parameters 



A = Vl + t'+ t> 



B = t'+ 
B* = t' . 



(70) 
(71) 



A* = Vl + f t'+ 
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